^^postgraduate  sew®* 

NAVAL  P0S1  VjtTTf     9?a*0 


6> 


-7 


AN  INVESTIGATION  OF  THE  PERFORMANCE 
AND  MIXING  PHENOMENA  ASSOCIATED  WITH  A  SUPERSONIC  EXHAUSTER 
INTERACTING  WITH  SUBSONIC  SECONDARY  FLOW 

by 


Philip  Alvah  Monroe 
Lieutenant  Commander,  United  States  Navy 
B.A. ,  Cornell  University,  1957 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 

AERONAUTICAL  ENGINEER 

from  the 


NAVAL  POSTGRADUATE  SCHOOL 
September  1S67 


ABSTRACT 

The  Transonic  Turbine  Test  Rig  at  the  Naval 
Postgraduate  School,  Monterey,  California,  uses  a  super- 
sonic jet  exhauster  to  lower  the  turbine  exit  static 
pressure  and  thus  increase  pressure  ratios  available  for 
turbine  testing.   Exhauster  operation  results  in  inter- 
action between  two  co-flowing  streams  of  air.   The  primary 
air  is  supersonic  and  the  secondary  air,  subsonic.   This 
thesis  is  an  experimental  study  of  the  mixing  process  that 
occurs  between  the  two  streams  in  a  cylindrical  mixing 
chamber.   In  addition,  a  theoretical  model  of  the  exhauster 
performance  based  on  poly tropic  efficiencies  and  conserva- 
tion of  momentum  was  programmed,  and  the  predicted  values 
were  compared  to  test  results. 


Thesis  by:   Philip  A.  Monroe  entitled  An  Investigation  of 
Performance  and  Mising  Phenomena  Associ .....  ..  Supers-.,  c 

Exhauster  Interactina  with  Subsonic  Secondary  Flow 
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1.   Introduction 

Exhausters  are  being  incorporated  into  the  design  of 
many  systems  because  of  the  simplicity  of  their  design  and 
the  absence  of  moving  parts.   Air  augmentation  of  rocket 
thrust,  torpedo  drag  reduction  innovations,  use  as  a  jet 
pump  in  aircraft  for  fuel  transfer  or  boundary  layer  con- 
trol are  significant  applications. 

The  Transonic  Turbine  Test  Rig  (TTTR)  installed  at  the 
Propulsion  Laboratory  of  the  Naval  Postgraduate  School  uses 
a  supersonic  exhauster  system  to  increase  the  range  of  tur- 
bine pressure  ratios  available  for  testing.   The  test  tur- 
bine is  enclosed  in  a  hood.   A  compressor  can  provide  a 
maximum  turbine  inlet  total  pressure  of  43  psia  giving  a 

maximum  turbine  pressure  ratio  (Pt     /P_    )  of  about  3 

inlet   exit 

when  exhausting  to  the  atmosphere.   Since  the  turbine  uses 

only  part  of  the  mass  flow  furnished  by  the  compressor,  the 

remainder  can  be  used  to  drive  the  exhauster  which  produces 

a  vacuum  in  the  hood  and  thereby  increases  the  turbine 

pressure  ratio. 

Very  little  information  is  available  on  the  mixing 

phenomena  encountered  with  supersonic  exhausters. 

Abramovich  provides  extensive  information  for  the  case  of 

an  incompressible,  free  jet  and  associated  turbulent  mixing! 
He  includes  test  work  accomplished,  and  theory  developed, 


Abramovich,  G.  N. ,  The  Theory  of  Turbulent  Jets 
(Cambridge,  MassachusettsT  The  MIT  Press,  Massachusetts 
Institute  of  Technology,  1963),  Part  I,  "Turbulent  Jets 
in  Incompressible  Fluid* " 


13 


by  Prandtl  and  Tollmien.   The  effect  of  compressibility  on 
turbulence  is  very  complicated  however,  and  simplifica- 
tions of  the  problem  must  be  introduced  to  develop  approx- 
imate solutions. 

This  thesis  is  a  study  of  the  supersonic  exhauster 
system  used  with  the  TTTR.   Objectives  of  the  investigation 
were  to  determine  the  hood  total  pressure  variation  with 
drive  nozzle  total  pressure  for  a  given  secondary  flow 
rate, and  to  study  the  mixing  and  turbulent  interchange  that 
occurs  between  the  two  co-flowing  streams  in  the  mixing 
pipe.   Main  topics  include:   traverse  probe  design,  experi- 
mental results  of  the  ejector  test  series,  a  computer  pre- 
diction program  of  exhauster  performance  based  on  poly- 
tropic  processes  and  conservation  of  momentum,  empirical 
determination  of  polytropic  efficiencies,  and  an  attempt  to 
determine  from  test  data  the  losses  due  to  mixing  and  tur- 
bulent interchange  between  the  two  streams  of  air. 

A  summary  of  the  general  properties  of  turbulent  jets, 
including  terminology  used  to  describe  the  different  mix- 
ing regions  that  occur  with  restricted  jet  operation  and 
co-flowing  streams,  is  given  in  Appendix  D.   The  terms 
ejector  and  exhauster  are  used  interchangeably  throughout 
this  thesis. 

The  author  is  most  appreciative  of  the  guidance,  coun- 
sel and  encouragement  provided  by  Professor  M.  H.  Vavra, 
the  laboratory  assistance  given  by  Mr.  J.  E.  Hammer,  and 
the  probe  fabrication  work  of  Mr.  T.  B.  Dunton. 
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2.   Installation. 

The  TTTR  installation  is  shown  schematically  in  Fig.  1, 
and  Fig.  2  is  a  photograph  of  the  turbine  hood  and  exhauster 
assembly.   High  pressure  air  is  supplied  by  an  Allis- 
Chalmers  VA  312  compressor  and  enters  Tank  1.   Air  to  the 
turbine  passes  from  Tank  1  through  a  calibrated  flow  mea- 
surement nozzle  and  enters  Tank  2.   The  amount  of  air  to  the 
turbine  is  controlled  by  a  manual  valve;  an  electrically 
operated  exhauster  inlet  valve  controls  the  drive  nozzle 
total  pressure  (PTN) .   More  detailed  descriptions  of  the 

TTTR  and  its  operating  performance  are  given  by  Eckert^  and 

3 
by  Commons . 

The  exhauster  section  of  the  TTTR  is  shown  in  greater 
detail  in  Fig.  3.   Important  features  are  the  arrangement 
of  twelve  mixing  tube  static  pressure  taps  (M2  -  M13) ,  the 
Kiel  probe  to  measure  drive  nozzle  total  temperature  (TTN) 
and  total  pressure,  and  the  secondary  air  static  pressure 
taps  just  prior  to  entry  into  the  mixing  pipe.   Secondary 
air  is  defined  as  the  air  that  enters  the  mixing  pipe  from 
the  turbine  hood.   To  obtain  the  mixing  pipe  static  pres- 
sures, four  taps  are  equally  spaced  circumferentially  at 


^Eckert,  R.  H. ,  Performance  Analysis  and  Initial  Tests 
of  a  Transonic  Turbine  Test  Rig  (Aeronautical  Engineer's 
Thesis,  U.  S.  Naval  Postgraduate  School,  Monterey, 
California,  May  1966),  Section  2,  "Installation." 

^Commons,  P.  M. ,  Instrumentation  of  the  Transonic 
Turbine  Test  Rig  to  Determine  the  Performance  of  Turbine 
Inlet  Guide  Vanes  Through  the  Application  of  the  Momentum 
and  Moment  of  Momentum  Equations  (Master's  Thesis,  U.  S. 
Naval  Postgraduate  School,  Monterey,  California,  1967) , 
Section  2,  "Installation." 
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each  location  and  commoned  to  obtain  an  average  reading. 
Secondary  air  static  pressures  are  measured  by  three  taps 
at  each  location,  equally  spaced  circuraferentially  and 
commoned. 

One  section  of  the  mixing  pipe  was  modified  by  cutting 
a  3  inch  by  1/4  inch  slit  at  3  locations  to  provide  ingress 
and  egress  of  a  traverse  probe.   The  slits  were  designed  to 
place  the  traverse  probe  total  pressure  opening  at  the  same 
distance  from  the  drive  nozzle  exit  as  the  mixing  pipe  wall 
static  pressure;  that  is,  in  the  plane  of  the  wall  static 
pressure  taps  for  a  given  axial  location.   Seals  were  also 
fabricated  to  prevent  leakages,  with  or  without  the  tra- 
verse probe  in  position. 

Secondary  air  total  temperature  (TTD)  is  measured  with 
a  modified  Kiel  probe.   Hood  total  pressure  (PTD)  is  mea- 
sured by  eight  static  pressure  ports  that  are  commoned  to 
obtain  an  average  value.   The  secondary  air  is  essentially 
that  which  passes  from  Tank  2  to  the  turbine.   It  is  aug- 
mented by  a  small  amount  that  enters  the  shaft  labyrinth 

seal  of  the  hood.   Dimensions  of  the  shaft  labyrinth  shown 

4 
in  Fig.  4  differ  from  those  given  by  Eckert  due  to  increas- 
ing shaft  clearance  from  0.005  inch  to  0.020  inch.  A 
thermocouple  is  placed  at  the  labyrinth  entrance  to  measure 
the  total  temperature  (TTC)  of  the  leakage  air. 

Drive  nozzle  I,  the  exhauster  drive  nozzle  designed  for 
the  TTTR,  is  shown  in  Fig.  5  with  detailed  measurements. 


4 

Eckert,  op.  cit. ,  Fig.  4,  p.  25. 
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Drive  nozzle  II,  a  smaller  nozzle  with  the  same  area  ratio 
(Atkroat/Aex^t) /  was  designed  and  built  to  investigate  the 
effects  of  different  primary  and  secondary  flow  areas, and 
smaller  primary  flow  rate, on  the  mixing  process  and  on  hood 
total  pressure.   Details  of  Drive  nozzle  II  are  shown  in 
Fig.  6.   Design  calculations  were  made  to  insure  a  uniformly 
decreasing  secondary  flow  area  for  the  new  nozzle. 

All  temperatures  are  measured  using  standard  Iron  - 
Constantan  or  Cromel  -  Alumel  thermocouples  with  an  ice  bath 
as  reference.   Differential  pressures  across  the  flow  mea- 
suring nozzle  are  measured  in  inches  of  water  on  a  mano- 
meter board  with  0.1  inch  gradations.   Mixing  tube  static 
pressures  are  measured  against  atmospheric  pressure  with  a 
Texas  Instruments  Fused  Quartz  Pressure  Gage.   All  other 
pressures  are  measured  in  inches  of  mercury  against  atmos- 
pheric pressure  as  reference  using  96  inch  long  manometer 
tubes  with  0.1  inch  gradations.   A  check  of  static  pressure 
readings  between  the  Texas  Instrument  gage  and  the  mercury 
manometers  produced  values  within  0.5  per  cent  of  each 
other. 

3.   Probe  Design. 

In  order  to  study  the  mixing  of  two  co-flowing  streams, 

velocity  traverses  were  desired  at  a  number  of  axial  loca- 

5 
tions.   Naviaux  attempted  to  obtain  this  information  but 


5Naviaux,  J.  C. ,  Transonic  Turbine  Test  Rig  Exhauster 
System  Tests  and  Tests  of  a  Reaction  Turbine  (Master's 
Thesis,  U.  S.  Naval  Postgraduate  School,  Monterey, 
California,  December  1966),  Section  4,  "Test  Results." 
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experienced  difficulty  in  obtaining  reliable  data.   During 
his  tests  the  probe  had  excessive  vibration,  the  static 
pressure  probe  used  had  the  pressure  ports  only  2  diameters 
away  from  the  probe  tip,  and  a  thermocouple  was  not  incor- 
porated/ hence,  the  differences  in  stagnation  temperatures 
between  the  two  streams  could  not  be  measured. 

For  the  ejector  tests,  the  four  probes  of  Fig.  7  were 
designed  and  built  to  obtain  traverse  information;  two  were 
pressure  probes  with  a  thermocouple  and  two  were  hot  wire 
probes.   Probe  I  was  the  only  successful  instrument  and  all 
traverse  information  was  obtained  with  it. 

Probe  I  is  shown  with  its  significant  dimensions  in 

g 

Fig.  8.   Criteria  presented  by  Gettelman  and  Krause  con- 
cerning static  port  location  for  maximum  accuracy  were 
strictly  adhered  to.   Best  results  are  obtained  when  the 
static  ports  are  4-6  sensing  element  diameters  away  from 
the  tip, and  greater  than  8  support  diameters  away  from  the 
support.   These  minimum  criteria  produce  static  pressure 
readings  within  5.0  per  cent  of  the  velocity  head  at 
Mach  0.9  and  3.5  per  cent  at  Mach  0.5.   Mr.  E.  Keener  of 
NASA  Ames  Laboratories  assured  the  author  that  the  error 
decreases  in  supersonic  conditions. 


6Gettelman,  C.  C.  and  Krause,  L.  N. ,  Considerations 
Entering  into  the  Selection  of  Probes  for  Pressure 
Measurement  in  Jet  Engines   (ISA  Proceedings  Vol.  7, 
Aeronautic  Instrumentation  -  Session  II),  pp.  134-137. 
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The  blockage  of  probe  I  is  4.3  per  cent  of  the  mixing 

7 
pipe  area.   Tests   show  that  a  blockage  effect  based  on  ratio 

of  support  diameter  to  duct  or  jet  diameter  can  induce  an 

error  as  large  as  8.0  per  cent  for  conditions  and  dimensions 

used  in  the  ejector  tests.   The  error,  however,  causes  a 

measurement  that  is  too  large  in  a  circular  duct  and  a 

measurement  that  is  too  small  in  a  jet.   This  difference  in 
direction  of  the  measurement  error  probably  leads  to  a 
partial  cancellation  of  the  blockage  error  at  mixing  stations 
near  the  drive  nozzle  for  the  exhauster  system.   The  exact 
effect  and  variations  with  axial  position  were  not  evaluated. 

The  thermocouple  of  probe  I  was  placed  in  the  support 
0.5  inch  above  the  pressure  element.   This;  position  was 
chosen  to  avoid  a  second  extended  sensing  arm  which  could 
induce  errors  due  to  probe  interactions „   The  thermocouple 
is  exposed  to  the  stream  by  a  0.160  inch  diameter  hole.   A 
0.040  inch  diameter  hole  in  the  after  side  provides  air 
circulation.   Total  temperature  measurements  at  mixing  pipe 
station  Ml  were  less  than  2.0  per  cent  different  than  those 
obtained  for  the  primary  and  secondary  flow  before  entry  to 
the  mixing  pipe. 

The  probe  vibration  problem  was  eliminated  by  instal- 
lation of  two  4  inch  x  6  inch  box  supports  contoured  to  fit 
the  outside  surface  of  the  mixing  pipe.   Each  box  was 


7 
Krause,  L.  N.  and  Gettelman,  C.  C„,  Effect  of  Inter- 
action Among  Probes,  Supports,  Ducts,  Walls  and  Jet 
Boundaries  on  Pressure  Measurements  in  Ducts  and  Jets   (ISA 
Proceedings  Vol.  7,  Aeronautic  Instrumentation  -  Session  II), 
Fig.  4,  p.  140. 
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fitted  with  a  3  inch  teflon  bushing  to  provide  probe  support 
over  a  long  distance.   This  eliminated  the  pivot  point  about 
which  a  probe  would  tend  to  oscillate.   The  assembly  is 
shown  in  detail  in  Fig.  9  and  pictured  in  Fig.  10.   The 
design  proved  very  successful. 

The  mixing  process  is  primarily  a  function  of  axial 
location  in  the  mixing  pipe.   Probe  I  has  a  3  inch  axial 
distance  between  the  position  where  the  total  pressure  is 
taken  and  the  thermocouple.   Because  significant  gradients 
in  total  temperature  exist,  particularly  near  the  drive 
nozzle,  probe  II  was  designed  to  obtain  the  temperature 
nearer  to  the  total  pressure  opening.   An  exposed  thermo- 
couple was  placed  in  a  position  1/2  inch  aft  of  the  static 
pressure  ports. 

Testing  of  probe  II  produced  expected  results  in  the 
subsonic  secondary  flow.   In  the  supersonic  region  of  the 
primary  flow,  the  test  data  were  irregular  both  for  static 
pressures  and  temperatures.   However,  lower  temperatures 
were  recorded  whereas  in  actuality  the  primary  flow  had  a 
higher  total  temperature  than  the  secondary  flow.   The  lower 
temperature  reading  occurred  because  an  exposed  thermocouple 

in  supersonic  flow  measures  an  adiabatic  wall  temperature. 

p 
Shapiro   introduces  the  energy  equation  to  describe  the 


8 
Shapiro,  A.  H. ,  The  Dynamics  and  Thermodynamics  of 

Compressible  Fluid  Flow  Volume  II   (New  York:  The  Ronald 

Press  Company,  1954) ,  pp.  1034-1036. 
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relationship  between  the  adiabatic  wall  or  recovery  tempera- 
ture  Taw,   the  static  temperature   Ts,   and  velocity  V 
that  must  exist  in  the  boundary  layer  over  such  an  exposed 
thermocouple . 


a-W 


R  is  the  recovery  factor  that  has  been  shown  in  many  experi- 

Q 

ments   to  be  equal  to  the  square  root  of  the  Prandtl  number 
or  0.85  in  laminar  flow.   Probe  II  could  be  used  in  super- 
sonic flows  only  if  the  velocity  were  known  and  actual  tem- 
perature, the  quantity  to  be  determined  from  measurement. 
This  temperature  was  needed  in  the  ejector  tests  to  deter- 
mine velocity  so  probe  II  was  not  applicable. 

An  attempt  was  made  to  verify  the  shape  of  the  pecul- 
iar velocity  profiles  obtained  at  axial  stations  Ml  -  M3 
by  using  a  hot  wire  anemometer.   The  results  of  such  a  test 
would  be  qualitative  because  of  the  temperature  difference 
between  the  primary  and  secondary  streams.   However,  the 
difference  in  total  temperature  never  exceeded  80°F.   With 
the  hot  wire  temperature  at  1200   F,  results  within  10  per 
cent  accuracy  could  still  be  obtained.   This  would  be  suffi- 
ciently accurate  to  confirm  the  velocity  profile  obtained 
with  probe  I. 

Probe  III  is  shown  in  Fig.  7.   It  is  a  hot  wire  mount- 
ed on  vertical  arms  attached  to  1/4  inch  tubing.   On  each 


9 
Ibid.,  pp.  1070  -  1072. 
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of  six  tests  the  wire  maintained  its  integrity  in  the  secon- 
dary flow  but  broke  as  soon  as  the  primary  jet  was  en- 
countered.  It  is  possible  that  the  hot  wire  broke  because 
of  vibrations  of  the  thin  vertical  arms  to  which  the  wires 
were  attached. 

Probe  IV,  shown  in  Fig.  7,  was  designed  such  that  the 
arms  to  which  the  wire  was  attached  were  stronger  and  ori- 
ented parallel  to  the  flow.   Both  changes  were  made  to 
dampen  arm  vibration.   The  probe  was  placed  in  the  center 
stream  before  the  exhauster  was  operated.   The  air  flows 
through  the  turbine  and  drive  nozzle  were  slowly  increased, 
but  the  wire  broke  at  velocities  far  below  those  experi- 
enced during  normal  exhauster  operating  conditions. 

4.   Description  of  Test  Runs. 

Tests  were  made  to  determine  the  hood  total  pressure 

variation  as  a  function  of  drive  nozzle  total  pressure  and 

secondary  flow  rate, to  establish  the  exhauster  performance, 

to  obtain  the  variation  of  mixing  pipe  static  pressure  as  a 

function  of  drive  nozzle  total  pressure  and  secondary  flow 

rate,  and  to  study  the  interaction  of  the  two  co-flowing 

streams  by  obtaining  velocity  profiles  at  different  axial 

positions  of  the  mixing  pipe.   Testing  was  accomplished 

with  two  different  drive  nozzles  in  order  to  compare  the 

effect  of  variation  of  the  primary  flow  rate  and  secondary 

flow  annulus  area  on  the  problem. 

A  total  of  34  runs  was  carried  out  between  February  and 
September  1967.   Exhauster  performance  and  mixing  pipe 
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static  pressure  rise  were  obtained  for  drive  nozzle  total 
pressure  ratios  (PTN/P  )  of  2.65,  2.45,  2.18,  1.90,  and 
1.63  and  for  secondary  flow  rates  of  4  and  2  lbm/sec. 

The  TTTR  was  set  in  the  following  manner.   Equation  1 
given  in  Sec.  5, entitled  "Data  Reduction",  is  the  equation 
used  to  determine  actual  flow  rate  through  the  flow  measur- 
ing nozzle.   Using  average  values  for   oc  ,   K,  and  Y, 
and  with  D2   known  for  the  TTTR  installation,  an  expression 
for  the  differential  pressure  1^  across  the  pressure  taps 
was  obtained. 


K  - 


—  IP1 

.2.18   V  r, 


w  =  desired  secondary  flow  rate 

p  =  static  pressure  at  upstream  tap 

T,  =  fluid  temperature  at  upstream  tap 

For  each  test  condition,  the  drive  nozzle  total  pressure  was 
set  by  hand.   After  allowing  time  for  stabilization  of  con- 
ditions, readings  of  the  calibrated  flow  measuring  nozzle 
total  temperature   T-,   and  pressure  p,   were  taken,  and  an 
hw  calculated  that  corresponds  to  the  desired  flow  rate. 
By  adjusting  the  turbine  inlet  valve,  this  value  of  h^ 
was  set.   Flow  rate  regulation  to  within  0.05  lbm/sec  of 
the  desired  amount  was  obtained. 

There  was  no  attempt  made  to  precisely  regulate  either 
primary  or  secondary  air  temperature  due  to  the  large  number 
of  variables  and  valve  settings  involved.   Turbine  testing 
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being  carried  out  concurrent  with  the  ejector  tests  created 
different  turbine  efficiencies  and  operating  conditions 
which  made  temperature  variations  difficult  to  regulate. 

The  analysis  of  early  test  results  showed  that  traverse 
probe  position  significantly  affected  hood  total  pressure 
and  secondary  flow  polytropic  efficiencies.   The  probe  was 
removed  during  final  tests  to  obtain  exhauster  performance 
and  mixing  pipe  static  pressure  rise. 

Velocity  traverses  were  initially  made  at  mixing  sta- 
tions Ml,  M2,  M3,  M4,  M6,  M7 ,  and  M8  of  Fig.  3  for  a  pres- 
sure ratio   PTN/Patrn  of  2. 18,  and  secondary  flow  rates  w. 
of  4,  3,  and  2  lbm/sec.   Data  points  were  taken  in  1/2  inch 
radial  increments  at  stations  M6  -  M8.   Due  to  large  gradi- 
ents in  total  temperature  and  pressure,  the  increment  was 
changed  to  1/4  inch  for  stations  M2  -  M4  and  to  1/10  inch 
at  Ml.   Additional  traverses  were  obtained  at  Ml  and  M4  for 
w   of  4   lbm/sec  and  PTN/P  .    of  2.65  and  1.63. 

It  was  not  possible  to  change  the  traverse  probe  posi- 
tion without  shutting  down  the  exhauster  system.   After 
restart  and  stabilization  of  the  system,  the  new  operating 
conditions  could  vary  in  total  temperature  up  to  3  degrees 
from  the  previous  run.   The  consequence  was  considered 
secondary  compared  to  the  effect  probe  position  had  on  ex- 
hauster operation. 

Testing  with  probe  II  in  run  17  produced  total  pres- 
sures different  from  those  obtained  with  probe  I  for  the 
same  operating  conditions.   Investigation  of  the  equipment 
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showed  a  leak  in  the  traverse  probe  total  pressure  tubing. 
The  effect  of  the  leak  became  noticeable  with  probe  II 
because  of  its  smaller  total  pressure  opening.   After  re- 
pairs, a  retest  of  data  points  indicated  the  leak  had 
existed  during  all  traverse  tests. 

The  planned  test  program  was  revised  and  retests  ini- 
tiated.  Final  tests  included  obtaining  the  hood  total 
pressure  and  mixing  pipe  static  pressure  rise  versus  drive 
nozzle  total  pressure,  for  secondary  flow  rates  of  4  and  2 
lbm/sec.   Velocity  traverses  were  taken  at  stations  Ml  -  M8 
(less  M5)  for  w.   of   4  and  PTN/Patm  of  2.18.   The  tests 
were  carried  out  for  both  drive  nozzles. 

Several  attempts  to  verify  the  velocity  profile  at 
station  Ml  with  a  hot  wire  anemometer  proved  futile.   The 
wire  could  not  withstand  the  impact  of  the  primary  stream 
and  broke  immediately  upon  insertion  in  every  test. 

5.   Data  Reduction. 

Secondary  Flow 

Secondary  flow  rates  were  measured  by  a  flow  nozzle 
between  tanks  1  and  2.   The  basic  formula  for  flow  rate  in 
lbm/sec  was  given  by  Eckert   to  be 


w     =    O.lliT  D^oc  KYJV2*  (1) 

r  i 


Eckert,  op_.  cit.  ,  Eq.  145,  p.  98, 


25 


D2  =  nozzle  throat  diameter  (in.) 

OC  =  coefficient  of  thermal  expansion  based  on 
2 
Dt.    (dimensionless) 

K  =  discharge  coefficient  (dimensionless) 

Y,  =  expansion  coefficient  (dimensionless) 

P1  =  static  pressure  at  upstream  tap  (in.  Hg.) 

T,  =  fluid  temperature  at  upstream  tap  (°R) 

1^  =  differential  pressure  across  taps  (in.  H2O) 

Values  for  OC     were  determined  as  outlined  by  Eckert. 
The  expansion  factor  Y,   was  changed  to  that  given  by  the 

ASME  Power  Test  Codes.  2   Discharge  coefficient   K  was 

13 
found  by  Naviaux   to  be  a  function  of  Reynolds  number. 

Naviaux  plotted  nozzle  flow  coefficient  versus  Reynolds 
number  and  by  curve  fitting  obtained  an  eighth  order  poly- 
nomial for   K,  with  Reynolds  number  the  independent  varia.- ". 
ble.   Computation  of  the  polynomial  showed  that  even  the 

highest  order  terms  produce  large  corrections  for   K.   For 

5 
Re  of  11.31  X  10   the  seventh  order  correction  was  +1168.857 


i:LIbid.  ,  pp.  98-99. 

12 

Flow  Measurement,  Supplement  to  ASME  Power  Test 

Codes,  New  York:  ASME,  1959,  Chap.  4,  Part  5,  p.  74. 

13 

Naviaux,  op.  cit.,  p.  35. 
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and  the  eighth  order  terra  -96.291.   Because  values  of   K 
were  determined  by  Naviaux  to  be  in  the  small  range  from 
1.00  to  1.05,  the  polynomial  with  such  large  corrections 
was  considered  inappropriate,  and   K  was  varied  in  steps 
of  0.01  as  a  function  of  Re   using  Naviaux' s  graphical 


results  as  a  reference 


14 


The  method  used  for  estimating  the  leak  rate  through 

15 
the  shaft  labyrinth  was  developed  by  Egli   and  revised  by 

16 
Eckert   to  apply  to  the  TTTR  installation. 
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OC   =  discharge  coefficient  for  a  single 
throttling  (dimensionless) 

^   =0.76  from  Egli  for  present  installation 

A   =  cross  sectional  area  of  seal  passage  (in.  ) 


14 


Ibid. ,  Fig.  19,  p.  68. 
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Egli,  A.  "The  Leakage  of  Steam  Through  Labyrinth 

Seals,"  Transactions  of  ASME,  V.  57,  1935,  pp.  115  -  122. 


16 


Eckert,  op.  cit. ,  pp.  103  -  105. 
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Pa.     =  total  pressure  at  seal  entrance  (psia) 

p   =  static  pressure  at  seal  exit  (psia) 

Tj.   =  total  temperature  at  seal  entrance  (°R) 

n   =  number  of  throttlings 
The  calculation  of  secondary  flow  rate  was  accomplished 
with  computer  program  SECFLO  which  is  shown  in  Appendix  C. 

Velocity 

The  velocity  profile  data  reduction  was  accomplished 
by  the  computer  program  VELOCITY.   Total  temperature,  total 
pressure  and  static  pressure  were  measured.   The  pressure 

ratio  was  tested  for  supersonic  or  subsonic  conditions. 

17 
For  supersonic  flow,  the  Rayleigh  Pitot  tube  formula   was 

used. 


pv 

c-v  kt 

% 
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For  subsonic  conditions,  the  standard  compressible  flow 
relation  was  used. 


f    -  ( I  +  ^  M 
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Shapiro,  A.  H. ,  The  Dynamics  and  Thermodynamics  of 

Compressible  Fluid  Flow  Volume' I   (New  York:  The  Ronald 

Press  Company,  1954) ,  p.  118. 
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In  either  case,  Mach  number  was  solved  by  iteration  and 
converted  to  velocity  with  the  measured  total  temperature, 
using: 


,2 


tt  =  ceT  +&      =>  V  =  y^Tcpfro-T) 


(7) 


The  computer  program  used  for  the  reduction  is  given  in 
Appendix  C. 

6.   Exhauster  Performance  Prediction. 

An  analysis  of  the  exhauster  performance  to  predict 
the  hood  total  pressure,  for  specified  drive  nozzle  total 
temperature  and  total  pressure,  and  secondary  air  mass  flow 
rate  and  total  temperature,  has  been  carried  out  by  Vavra-^-^ 

in  the  TTTR  preliminary  design  calculations,  and  used  by 

19 
Naviaux.     However,  both  analyses  assume  that  the  drive 

nozzle  is  operating  at  design  conditions,  and  the  static 

pressure  at  the  entrance  to  the  mixing  section  is  constant 

in  radial  direction.   Considering  the  many  variables  in  the 


18 

Vavra,  M.  H. ,  Theoretical  Evaluation  of  Exhauster 

System  of  Transonic  Turbine  Test  Rig   (TN65T,  U.S.  Naval 

Postgraduate  School,  Monterey,  California,  1965) ,  pp.  16 

21. 

19 

Naviaux,  op.  cit. ,  pp.  22  -27. 
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physical  set-up  of  the  TTTR,  such  as  turbine  efficiency, 
RPM,  secondary  air  mass  flow  rate  and  total  temperature, 
and  the  range  of  drive  nozzle  total  pressures  available,  it 
seems  unlikely  that  the  nozzle  will  operate  always  at  de- 
sign conditions.   Test  data  plotted  in  Fig.  29  indicate  a 
static  pressure  gradient  at  station  Ml,  located  just  9/16 
inch  from  the  nozzle  exit.   Therefore,  the  possibility  of 
a  discontinuity  in  static  pressure  between  the  nozzle  exit 
(area  Al)  and  the  secondary  flow  annulus  (area  AIT)  at  the 

start  of  the  mixing  section  is  contained  in  the  present 

20 
development.   This  hypothesis  is  supported  by  Jaumotte 

who  indicated  that  the  classic  hypothesis  that  assumes 
uniformity  of  the  static  pressure  in  the  initial  section  of 
mixing  is  wrong.   In  the  cited  reference  it  is  stated  also 
that  the  laws  governing  the  pressures  in  the  two  fluids  and 
the  factors  which  influence  them  are  not  yet  well  under- 
stood. 

The  complete  revised  analysis  of  exhauster  performance 
is  included  as  Appendix  A.   Appendix  B  gives  the  computer 
program  used  to  calculate  prediction  maps  of  exhauster  per- 
formance over  its  operating  range.   The  pertinent  steps  of 
the  analysis  are  outlined  below  with  the  notations  corres- 
ponding to  those  of  Fig.  3. 


20 

Engle,  M.  0. ,  Some  Problems  in  the  Design  and 

Operation  of  Jet  Ejectors,  Proceedings  of  the  Institution 
of  Mechanical  Engineers  Steam  Plant  Group   (London,  England, 
No.  133,  1963),  p.  350. 
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Polytropic  process  relations  as  well  as  non -dimensional 

21 

flow  functions  are  used.   These  are  taken  from  Vavra.     For 

an  expansion  process 


wV^ 


/U.   it 


I    '' 


where 


n  = 


n   =  polytropic  exponent 

dT 
T2p  =  polytropic  efficiency    , 
f  aTis 

0   =  ratio  of  specific  heats 
w   =  flow  rate  (lbm/sec) 

Tx.   =  inlet  stagnation  temperature  (°R) 
^o 

A2   =  flow  area  at  station  of  interest  (in.'') 
Pi.   =  inlet  stagnation  pressure  (psia) 

P2   =  static  pressure  at  station  of  interest 
(psia) 

g   =  universal  gravitation  constant 
(lbm- ft/lb -sec2) 

Subscript  1  refers  to  inlet  conditions 
Subscript  2  refers  to  station  of  interest 


2. 
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Vavra,  M.  H. ,  Problems  of  Fluid  Mechanics  in  Radial 

Turbomachines  (Rhode-Saint-Genese,  Belgium:   Von  Karman 

Institute  for  Fluid  Dynamics,  1965) ,  VKI  Course  Note  55a, 

pp.  22  -  35. 
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For   a  compression  process 
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P^  =  static  pressure  at  1  (psia) 
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For  a  choked  nozzle 
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Initially,  polytropic  efficiencies  were  assumed  to  be 
constant  over  the  operating  range  of  the  exhauster.   The 
determination  or  estimation  of  the  values  used  is  described 
in  Sec.  8  entitled  "Determination  of  Polytropic  Exponents." 

Drive  Nozzle 


The  flow  functions  from  the  nozzle  entrance  to  the 
throat  (  $ n  )  ,  and  from  the  nozzle  entrance  to  the  exit 
(  §^  ) i  are 
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The  relationship 


(13) 
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Expressing  the  exit  flow  function  in  terms  of  a  presi- 
sure  ratio 


ex 


(15) 


With   ^     given  by  Eq.  14  and  PTN  an  initial  condi- 


ex 


tion,  PIN  may  be  solved  by  iteration.   With  PIN,  the 
quantities   TIN,  VlN,  and  w   may  then  be  obtained. 


Secondary  Flow 

For  the  expansion  process  from  the  hood  to  the  mixing 
pipe  annulus  area,  there  is 


B.    =  £ 


zy 


Jf-I  Llli 


f,T  \V" 


/>j  ♦ 


^      -    41 '  Pu 


(16) 


TO 


With  secondary  flow  rate  and  hood  total  temperature 
known,  a  value  of  PIT  is  assumed  and  PTD  calculated  by 
iteration  and  matching.   For  the  assumed  PIT  and  correspond- 
ing PTD,  the  values  of  TIT  and  VlT  can  be  calculated. 
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Mixing  Pipe 
With  the  values  calculated  for  the  drive  nozzle  and  for 
assumed  secondary  air  flow  conditions,  the  momentum  equation 
was  used  for  the  region  within  the  mixing  pipe  to  obtain 
conditions  at  the  diffuser  entrance.   The  basic  equation  is 

(&ii)yt  -  ftN  -±v„  =  (PirjAiT  +  (pin)AI  -(P2)A2  -  Ff       (17) 

Assumptions  include  uniform  velocity  profiles,  uniform  sta- 
tic pressure  distribution  over  each  area,  and  adiabatic 
flow.   The  term  Ff   represents  the  viscous  shear  force 
imposed  by  the  wall  on  the  moving  fluid.   Estimation  of  the 
magnitude  of  this  force  is  outlined  by  Naviaux.   In  the 
calculation,  the  resultant  pressure  drop  produced  by  Ff 
is  determined  by  a  single  iteration  of  Eq.  17. 

In  Eq.  17,  P2   and   V2   are  unknown.   With  continuity, 
the  pressure   P2   can  be  expressed  with  the  known  quanti- 
ties and  V2/  which  leaves  one  equation  to  be  solved  with 
the  unknown  quantity  V2.   With  V*2,  all  conditions  at  the 
diffuser  entrance  can  be  obtained. 

Diffuser 
For  the  compression  in  the  subsonic  diffuser 


mn-i=mm"*-im  ,i" 
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The  pressure   P3   is  obtained  by  iteration;  and  if  it  does 
not  match  atmospheric  conditions,  a  new  value  of  PIT  for  the 
secondary  air  must  be  assumed  and  the  calculation  repeated. 

7  o   Test  Results  and  Discussion. 

Hood  Total  Pressure 

Tests  of  the  exhauster  system  were  conducted  to  measure 
the  hood  total  pressure  obtained  over  the  operating  range  of 
the  drive  nozzle.   Runs  were  made  at  a  secondary  flow  rate 
of  4  and  2  lbm/sec  for  each  drive  nozzle. 

Two  attempts  to  obtain  the  information  were  discarded. 
Run  11  was  made  with  the  traverse  probe  at  station  M2  of 
the  mixing  pipe.   Comparison  with  other  test  data  showed 
that  hood  total  pressure  was  increased  by  as  much  as  1  inch 
of  mercury  with  the  traverse  probe  installed,  particularly 
when  the  probe  was  close  to  the  drive  nozzle.   The  increase 
in  hood  pressure  is  the  result  of  losses  due  to  shock  waves 
that  attach  to  the  probe,  the  drag  of  the  probe  in  the 
stream,  and  the  wake  formed  behind  the  probe.   Run  29  with 
drive  nozzle  II  was  not  used  due  to  the  discovery  that  a 
large  seal  in  the  hood  had  not  been  properly  installed. 
Run  18  and  22  for  drive  nozzle  I  and  run  32  for  drive  nozzle 
II  were  successful.   The  results  are  plotted  in  Fig.  11. 

The  graph  shows  that  for  the  same  secondary  flow  rate, 
drive  nozzle  I  creates  a  lower  hood  total  pressure  than 
nozzle  II.   This  occurs  primarily  because  nozzle  I  has  a 
throat  area  of  10.42  in.   and  nozzle  II  a  throat  area  of 
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8o29  in.2  so  that  with  choked  conditions,  about  1.51  lbm/sec 
greater  flow  rate  of  primary  air  exists  with  nozzle  I. 
Hence,  its  driving  momentum  flow  is  larger. 

For  drive  nozzle  I,  a  lower  hood  pressure  resulted  for 
a  smaller  secondary  flow  rate  at  a  given  drive  nozzle  pres- 
sure ratio.  In  a  relative  sense,  the  effect  is  the  same  as 
described  above.  For  the  same  drive  nozzle  pressure,  which 
establishes  the  mass  flow  rate  and  momentum  flow  of  the 
primary  nozzle,  it  is  possible  to  produce  lower  pressures 
in  the  hood  for  smaller  secondary  flow  rates. 

Mixing  Pipe  Static  Pressure  Rise 
Mixing  pipe  static  pressures  were  recorded  for  each 
pressure  ratio   PTN/patm  anc^  fl°w  rate.   Graphs  showing 
the  static  pressure  variation  as  a  function  of  distance 
from  the  drive  nozzle  exit  are  given  in  Figs.  12,  13,  and 
14.   In  general,  there  occurs  a  rapid  increase  in  the  sta- 
tic pressure  along  the  initial  region  of  the  mixing  pipe. 
In  the  main  mixing  zone  where  the  primary  and  secondary 
core  regions  have  disappeared,  there  is  a  slight  but  steady 
decrease  in  static  pressure  in  each  case,  except  for 
PTN/Patm  =  2.65  in  Fig.  12.   The  exception  could  be  due  to 
the  fact  that  incomplete  mixing  has  occurred  in  the  mixing 
pipe  for  the  extreme  condition  tested.   This  point  is  the 
greatest  drive  nozzle  pressure  ratio  used  and  lowest  hood 
total  pressure  produced  during  testing.   As  expected,  the 
static  pressure  rises  in  the  diffuser  to  the  atmospheric 
pressure  at  the  exit. 
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No  attempt  was  made  to  predict  analytically  the  mixing 
pipe  static  pressure  as  a  function  of  axial  location.   Due 
to  the  extremely  complex  nature  of  the  flow  which  includes 
mixing,  turbulence,  and  shocks  in  this  region,  the  momentum 
equation  used  in  the  exhauster  prediction  analysis  is 
applied  only  for  the  whole  mixing  pipe  to  obtain  average 
conditions  ahead  of  the  diffuser.   The  rapid  rise  of  static 
pressure  in  the  initial  region  is  due  to  the  slowing  down 
of  the  supersonic  stream  as  it  mixes  with  the  subsonic 
secondary  air.   Also  present  is  a  series  of  expansions  and 
shock  waves  to  equalize  the  static  pressure  between  the  two 
flows.   A  possible  extension  of  this  thesis  would  be  to 
devise  a  mixing  chamber  for  Schlieren  photographic  study 
of  the  mixing  phenomena  in  compressible  flow. 

The  slight  but  steady  decrease  in  static  pressure  in 
the  main  mixing  zone  indicates  the  loss  of  available  energy 
due  to  losses  caused  by  the  shear  forces  in  the  flow  and 
the  continued  turbulent  interchange  between  particles. 

Figures  13  and  14  show  that  at  L/D  of  about  4,  the 
wall  static  pressure  is  higher  for  a  drive  nozzle  total 
pressure  ratio  of  2.45  than  for  2.18.   This  condition  re- 
peated itself  on  4  different  occasions,  but  each  time  the 
manometer  board  was  very  unstable  with  pressure  variations 
at  stations  Ml  —  M3  up  to  plus  or  minus  3  in.  Hg.  observed 
out  of  6  -  8  in.  Hg.  vacuum.   The  unsteady  condition  was 
not  reflected  in  the  hood  total  pressure.   It  appeared  to 
be  the  result  of  a  shock  moving  fore  and  aft  in  the  initial 
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region  of  the  mixing  pipe.   A  very  definite  sound  variation 
was  heard  while  operating  at  this  condition  which  would 
tend  to  support  the  hypothesis  of  an  unstable  shock  move- 
ment. 

Velocity  Profiles 
Velocity  profiles  obtained  with  probe  I  for  drive 
nozzle  I  are  shown  in  Figs.  15  through  19,  and  those  for 
drive  nozzle  II  in  Fig.  20.   There  occur  large  gradients 
between  the  two  co -flowing  streams  at  station  Ml  and  con- 
siderable non-uniformities  of  the  velocity  profiles.   Both 
of  these  effects  decrease  with  increase  in  axial  distance 
from  the  drive  nozzle. 

The  non-uniformity  of  velocity  profile  had  been  ex- 
po 
pected  because  of  Naviaux's  results   but  not  to  the  degree 

observed  in  the  ejector  tests.   However,  Naviaux  did  not 
traverse  closer  than  4  inches  to  the  drive  nozzle,  whereas 
at  station  Ml  the  traverse  total  pressure  opening  is  only 
9/16  inch  from  the  drive  nozzle  exit  plane. 

The  unusual  profiles  are  probably  caused  by  the  shock 
structure  that  actually  exists  in  the  flow.   A  full  discus- 
sion with  sketches  of  the  possible  pattern  is  included  in 
Sec.  10.   Essentially,  the  presence  of  internal  oblique  and 
normal  shocks  causes  sharp  gradients  in  pressure  and  veloc- 
ity in  the  flow  with  accompanied  differences  in  flow  direc- 
tions.  The  large  peak  in  static  pressure  shown  in  Fig.  29 


22 

Naviaux,  op.  cit. ,  Fig.  3,  p.  52. 
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for  station  Ml  is  probably  due  to  a  large  radial  velocity 

component  in  the  mixing  zone  at  this  radial  position. 

23 
Abramovich  cites  experiments  by  Zhestkov   for  a  free 

jet  operating  at  off-design  conditions.   The  results  are 

similar  to  the  ejector  tests  with  regard  to  the  decrease 

in  velocity  in  the  center  of  the  primary  jet, compared  to 

that  in  the  jet's  boundaries.   By  comparing  the  velocity 

profiles  at  stations  Ml,  M2,  and  M3  it  is  apparent  that  the 

secondary  air  increases  in  velocity  and  the  secondary  core 

region  decreases  in  size.   The  primary  stream  velocity 

decreases  and  the  primary  core  area  decreases  as  the  mixing 

zone  increases  in  size. 

Interesting  results  were  obtained  by  non-< 

dimensional izing  the  velocity  profiles  obtained  in  the 

main  mixing  region  of  the  mixing  pipe, in  the  manner  given 

24 
by  Abramovich   for  compressible  flows.   The  main  mixing 

region  is  the  part  of  the  mixing  pipe  where  the  "core" 

flows  of  primary  and  secondary  air  have  disappeared. 

Appendix  D  describes  the  regions  more  fully.   Abramovich 

uses  the  ordinate  V/Vax  where  V  is  the  velocity  at  a  given 

location,  Vax  is  the  axial  velocity.   For  the  absissa, 

Y 

2.2T  yc 


23 

Abramovich,  op.  cit. ,  pp.  335  -  338. 

24Ibid. ,  pp.  275  -  276. 
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where  y   and  yc   represent  the  coordinate  of  the  point 
which  corresponds  to  an  arbitrary  value  of  the  diinens ion- 
less  velocity, and  the  coordinate  where  the  velocity  is  1/2 
the  axial  velocity,  respectfully.   Figures  21  and  22  graph 
the  result.   In  both  cases  a  single  profile  is  obtained. 

Station  M4  is  not  inside  the  main  mixing  region  with 
drive  nozzle  I, but  it  is  with  drive  nozzle  II.   Figure  23 
shows  that  the  non-dimensional  profile  at  station  M4  for 
nozzle  I  deviates  from  the  single  profile  obtained  by  non- 
dimensionalizing  the  velocity  profiles  of  stations  M6  -  M8. 
From  these  results,  it  appears  that  a  sufficient  condition 
to  determine  the  start  of  the  main  mixing  region  in  a 
restricted  jet  is  to  determine  where  the  velocity  profile 
may  be  non-dimensional ized  to  a  single  profile.   The 
following  paragraphs  discuss  the  mixing  regions  in  greater 
detail. 

Mixing  Zone 
The  shape  of  the  mixing  zone  determined  from  the  velo- 
city profiles  is  shown  in  Figs.  24  and  25  for  nozzle  I  and 
nozzle  II,  respectively.   This  is  at  best  an  approximation 
because  of  the  limited  number  of  traverse  positions  in  the 
initial  mixing  region.   The  mixing  zone  appears  to  touch 
the  mixing  pipe  wall  at  an  L/D  of  3.78  for  nozzle  I, and 
4.71  for  nozzle  II.   D   is  the  nozzle  exit  diameter. 
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25 
Victorin   carried  out  a  series  of  experiments  using 

water  as  the  driven  and  driving  fluid  and  determined  that 

the  velocity  distribution  in  the  mixing  chamber  behaved  in 

the  same  fashion  as  a  circular  wake  with  the  width  propor- 

1/3 
txonal  to  x    ,  where  x   is  the  axial  position.   The 

assumption  is  made  that  the  primary  fluid  emanates  at  a 
point  source  at  x  =  0   located  at  a  position  upstream  of 
the  nozzle  so  that  the  velocity  distribution  would  have 
spread  to  match  the  nozzle  exit  dimensions.   Calculations 
of  the  theoretical  distance  x  '  ,  to  predict  where  the 
mixing  zone  meets  the  mixing  pipe  wall,  gave  an  L/D  of  8.27 
for  nozzle  I,  and  9.71  for  nozzle  II.   These  conditions  are 
shown  in  Fig.  26. 

It  appears  from  the  traverse  information  that  the  out- 
ward spread  of  the  mixing  zone  occurs  at  a  faster  rate  than 
the  incompressible  approximation.   From  the  test  data  ob- 
tained, a  closer  approximation  for  the  spread  of  a  jet  in 

2/5 
compressible  flow  would  be  y  =  x   .   The  relation  predicts 

that  mixing  zone  touches  the  pipe  wall  at  an  L/D  of  3.87 
for  nozzle  I, and  4.86  for  nozzle  II. 

From  the  traverse  data  from  puns  4  -  17  it  can  be  con- 
cluded that  the  shape  of  the  mixing  zone  is  independent  of 
secondary  air  mass  flow  rate.   This  conclusion  is  valid 
even  though  a  leak  existed  in  the  traverse  probe  total 


^Victorin,  k.  ,  Investigation  of  Turbulent  Mixing 
Processes,  NACA  TM  10  96,  1946. 
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pressure  line  during  these  runs.   The  leak  did  not  affect 
the  radial  position  of  the  steep  total  pressure  gradient 
observed  in  test  data.   The  radial  position  where  the  steep 
gradient  in  total  pressure  occurred  was  invariant  as  the 
secondary  flow  rate  was  set  at  4,  3,  and  2  lbm/sec. 

8.   Determination  of  Polytropic  Exponents. 

General 
The  theoretical  model  which  has  been  assumed  for  the 
prediction  of  exhauster  performance  uses  flow  functions 
with  polytropic  exponents  to  describe  the  flow  in  four 
regions  of  the  exhauster;  namely,  nozzle  entrance  to  nozzle 
throat  (nt) ,  nozzle  entrance  to  nozzle  exit  (n, ) ,  secondary 
air  from  the  hood  to  the  annulus  at  the  start  of  the  mixing 
pipe  (n^) ,  and  the  compression  process  in  the  diffuser 
(nc) .   The  relationship  between  polytropic  efficiencies 
(  Y?      )  and  polytropic  exponents  (n)  ,  for  an  expansion  pro- 


r±   = 


*-nj*->) 


b   ~  U-'  A  r)e  / 


and  for  a  compression  process 
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26 
Eckert   assumed  polytropic  efficiencies  to  be: 

TlT   =  0.975,    72  -  0-925,  T^    =  0.90,  and   T)     =  0.70, 


Secondary  Flow 
Tests  were  conducted  to  determine   n^   and  T£  .   for 
the  secondary  flow.   Secondary  flow  rate  (w.)  was  maintained 
at  4  lbm/sec,  and  drive  nozzle  total  pressure  ratio 
PTN/Patm  was  set  at  2.65,  2.45,  2.18,  1.90,  and  1.63.   The 
flow  function  for  this  expansion  is 


ttJBpTJT  -  $ 


(19) 


Quantities  on  the  left  were  either  measured  or  are  con- 
stants.  The  pressures  on  the  right  were  measured,   n^  was 
determined  by  iteration  in  computer  program  PDLY.   Results 
of  the  tests  are  given  in  Table  I.   For  4  lbm/sec  secondary 
flow  rate, the  average  values  are   n^  =  1.32  and  T)       =  0.85 
for  the  supersonic  operating  range  of  the  exhauster.   Con- 
ditions during  the  test  were  stable. 

Initial  testing  to  determine   n^   for  secondary  flow 
rate  of  2  lbm/sec  was  not  immediately  successful.   In  two 
separate  tests,  results  showed  polytropic  exponents  from 
1.08  up  to  1.41, and  one  value  of  1.61  which  is  impossible 


26 

Eckert,  op.  cit. ,  p.  93. 


43 


since   nd   cannot  be  larger  than   X     =1=41.   During  the 
test,  annulus  static  pressure  was  unsteady  and  an  average 
reading  was  used  for  data  reduction.   Figure  13  indicates 
that  an  unusual  condition  exists  in  the  mixing  pipe  between 
pressure  ratios  of  2,45  and  2„18„   The  curves  showing  mix- 
ing pipe  static  pressures  cross  each  other  for  the  given 
conditions.   This  phenomenon  repeated  itself  on  four  sepa- 
rate occasions. 

In  order  to  determine  the  effect  of  a  small  error  in 
the  annulus  static  pressure  readings,  the  computer  program 
POLY  was  used.   Values  for  hood  total  temperature,  total 
pressure,  and  secondary  flow  rate  were  held  constant  and 
values  of  annulus  static  pressure  varied  in  increments  of 
0.1  inch  mercury  or  0.05  psia.   The  result  is  given  in 
Table  II.   It  is  evident  that  annulus  static  pressure  is 
critical  to  the  calculation  of  nd   since  a  change  of 
JL.  0.1  psia  from  the  test  value  changes   n^  between  1.39 
and  1.24.   In  the  final  test  it  was  attempted  to  take  all 
readings  simultaneously.   Results  of  this  test  are  con- 
sidered accurate  and  are  given  in  Table  III.   They  show 
that  nd   is  about  1.28,  with  T^  d  equal  to  0.84.   The 
tests  showed  that  the  poly tropic  efficiency  was  independent 
of  drive  nozzle  total  pressure  but  dependent  on  secondary 
flow  rate. 

Dif fuser 
Diffuser  polytropic  efficiency  was  determined  empiri- 
cally from  the  static  pressure  rise  between  the  mixing  pipe 
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station  M9  to  the  diffuser  exit  for  all  operating  conditions 
tested.   Figure  27  is  a  graph  of  the  diffuser  static  pres- 
sure rise  for  secondary  flow  rates  of  4  and  2  lbm/sec  for 
different  pressure  ratios  of  drive  nozzle  total  pressure  and 

atmospheric  pressure. 

Initially,  the  exhauster  performance  prediction  program 
was  used  with  the  actual  exhauster  operating  conditions  as 
inputs  and  with  a  diffuser  efficiency   7^   =  0.7  as  assumed 
by  Eckert.   The  predicted  static  pressure  rise  from  M7  to 
atmospheric  conditions  based  on  Y[       of  0.70  is  plotted  also 

in  Fig.  27  for  4  and  2  lbm/sec  secondary  flow,  showing  that 
the  diffuser  efficiency  is  larger  than  0.7.   A  diffuser 

efficiency  of  0.90  was  assumed,  and  the  results  of  the  cal- 
culated static  pressure  rise  in  the  diffuser  is  shown  in 
Fig.  28.   Agreement  between  calculated  and  actual  diffuser 
static  pressure  rise  is  good,  but  the  graph  shows  that  there 
is  a  small  variation  of  diffuser  efficiency  with  drive  noz- 
zle total  pressure  ratio.   A  diffuser  efficiency  of  0.90 
very  accurately  predicts  the  performance  at  the  lower  drive 
nozzle  total  pressure  ratios  up  to  about  2.1  but  under- 
estimates the  performance  at  the  higher  pressure  ratios. 
The  prediction  program  was  altered  to  vary  T[c      with  drive 
nozzle  total  pressure.   Values  used  for  diffuser  efficiency 
at  different  drive  nozzle  total  pressure  ratios  are  shown 
in  Table  IV.   The  results  of  calculated  pressure  rise  in 
the  diffuser  compared  to  test  data  are  in  Fig.  28.   The 
curves  show  good  agreement  between  calculated  and  test 
results. 
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It  was  noted  also  that  the  diffuser  efficiency  varied 
slightly  with  secondary  flow  rate.   This  effect  is  small 
and  need  not  be  taken  into  account  in  the  simplified  model 
used  for  the  prediction  program. 

9.   Determination  of  Average  Shear  and  Friction  Coefficient. 

Mathematical  solutions  do  not  exist  that  exactly  pre- 
dict the  conditions  in  a  flow  as  complicated  as  that  in  the 
mixing  tube  of  the  exhauster.   However,  it  was  attempted  to 
establish  a  value  for  the  average  shear  stresses  in  the 
flow  by  an  approximate  method.   Such  an  approximation  was 
obtained  by  not  attempting  to  satisfy  the  hydrodynamic  dif- 
ferential equations  for  each  particle  but  by  integrating 
the  profile  distributions  obtained  in  laboratory  tests  and 
hence  solving  the  momentum  equation  for  the  flow  between 
two  stations.   The  momentum  equation 

U.U,    +    (IJJA,    -      fVicU^    +     J  TJJi^+TA,  (20) 

©  J©  J®  is> 

where:     1  represents  conditions  at  the  upstream  station 

2  represents  conditions  at  the  downstream  station 

A  =  surface  area  of  the  mixing  pipe  between 
stations 

can  be  used,  but  in  turbulent  flow  with  large  velocity  gra- 
dients T   consists  of  the  "true  shear"  due  to  the  visco- 
sity of  the  fluid,  and  the  "apparent  shear"  due  to  the  losses 
produced  by  the  turbulent  interchange  of  fluid  particles. 
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27 
Prandtl  states    that  except  in  layers  directly  next  to  the 

wall,  the  true  shear  is  small  enough  to  be  neglected  in  com- 
parison to  the  apparent  shear.   The  approach  assumes  axi- 
symmetric  flow,  a  perfect  gas,  axial  velocities  everywhere 
and  steady  flow,  with  gravitational  effects  being  ignored. 
Equation  20  is  used  for  the  profiles  of  velocity,  pressure 
(total  and  static) ,  and  temperature  measured  at  stations 
Ml,  M4,  and  M8. 

With  the  pressures  in  psi  and   r   in  inches,  there  are 


<AA  r  2.TT  rAr 


Equation   20  becomes 


BIT       f.V, 
•3 


Numerical  integration  of  large  scale  reproductions  of  Figs 
29  and  30  gave  at  station  Ml 


2IE  I  *V,2  i    . 


frJir  =  34Z.5"  +  473". O   -  817.5"  Ids 
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Prandtl,  L. ,  The  Mechanics  of  Viscous  Fluids, 

Aerodynamic  Theory   (California:  Durand  Reprinting 

Committee,  California  Institute  of  Technology,  January 
1943),  Vol.  Ill  G,  p.  135. 
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at  M4 

3*     J     ^ 

and   at  M8 


fr<kr    =    33HM     t  SZf.f    s     849.^     IV*. 


® 
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The  results,  which  show  an  increase  in  momentum  from 
station  Ml  to  M4,  must  be  in  error  although  use  of  the 
momentum  equation  to  model  the  mixing  pipe  performance  and 
the  mathematical  integration  can  be  considered  correct. 

The  error  could  have  been  introduced  by  the  blockage  effect 
of  the  traverse  probe  and  the  variation  of  the  overall 

exhauster  performance  with  probe  position.   To  illustrate 
this  condition,  Table  V  shows  the  values  of  the  hood  total 
pressure  for  the  same  exhauster  operating  conditions,  but 
with  the  traverse  probe  at  different  locations.   A  varia- 
tion of  1.65  in.  Hg.  hood  total  pressure  occurs  between 
the  conditions  where  the  traverse  probe  is  at  station  Ml, 
and  the  probe  is  completely  removed.   Another  possible 
reason  for  the  error  is  that  it  was  not  possible  to  produce 
exactly  the  same  settings  of  pressures  and  temperatures  for 
each  traverse  since,  to  change  traverse  probe  locations,  it 

was  necessary  to  shut  down  the  turbine  and  exhauster .   Re- 
establishment  of  all  previous  conditions  was  extremely 
difficult.   Analysis  of  the  data  indicates,  however,  that 
this  effect  should  be  secondary. 
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The  calculation  was  continued  between  station  M4  and 
M8  using  the  integrated  values  obtained. 

fj    =     869.1  -  787.3    -    82.6     lb,. 
f\s  =     6.  *Y6     -FT* 


From  the  test  data  and  results  obtained,  the  turbulent 

friction  coefficient   f  was  also  calculated.   Using 

28 
Vavra,    the  pressure  drop   A  p  due  to  friction  may  be 

approximated  by 

AT3  =  t  "pi f"Psf ) 


where     f  =  turbulent  friction  coefficient  (dimension- 
less) 

L  =  length  of  pipe  (in.) 

D  =  diameter  of  pipe  (in.) 

P   =  average  density  between  stations  (slugs/ft3) 

V  =  average  velocity  between  stations  (ft/sec) 

with 

p  -   rt  §     "■  too-  y       (slugs /ft  ; 

4 .  m.l  fn 


28 

Vavra,  M.  H. ,  Theoretical  Evaluation  of  Exhauster 

System  of  Transonic  Turbine  Test  Rig   (TN65T,  U.  S.  Naval 

Postgraduate  School,  Monterey,  California,  1965) ,  p.  26. 
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The  total  frictional  force  may  be  expressed  in  terms  of  a 
pressure  drop 

(J  =  Tf  77T  *?     =0-303  *f    (lbs.) 
With  Ff  determined,  f  becomes 


*  =  58ir  lTv*  (21> 


Average   values    used   to   calculate      f     were:      T  =   530   R, 
p  =   12.4  psia,    and  V  =  700    ft/sec. 


f  =   0.106 

This  is  one  order  of  magnitude  larger  than  the  friction 

29 
factor  0.017  predicted  by  the  Moody  diagram   for  commer- 
cial steel  pipes.   The  value  should  have  been  closer  even 

with  velocities  of  700  ft/sec.   Shapiro  cites  numerous 

30 
tests    that  confirm  the  use  of  the  incompressible  friction 

formula  for  conditions  up  to  Mach  =  1.0. 

The  calculations  were  repeated  between  Ml  and  M8. 

With  Ff  =  29.7  lbs.,  L  =  86  in.,  p  =  11.3  psia,  T  =  500°R, 

and  V  =  750  ft/sec,  Eq.  21  gives 

f  =  0.024 
This  value  is  still  large  but  has  the  correct  order  of 
magnitude  for  the  friction  coefficient. 


29 

Shapiro,  op.  cit. ,  Fig.  6.15,  p.  184 

30Ibid. ,  p.  1071. 
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The  results  of  the  calculations  for  friction  coeffi- 
cient suggest  that  the  traverse  data  at  station  M4  are  not 
accurate.   Better  results  were  obtained  between  Ml  and  M8 
than  M4  and  M8.   No  particular  quantity  in  the  M4  data  was 
considered  unreasonable.   The  conclusion  is  that  the  cal- 
culation of  an  average  shear  force  and  friction  coefficient 
for  the  mixing  pipe, using  the  ejector  test  data, produced 
unreliable  results. 

10.   Evaluation  and  Discussion  of  Predicted  Results. 

Actual  test  conditions  were  input  to  the  exhauster 
prediction  program  to  determine  a  predicted  hood  total 
pressure.   Three  predictions  were  made:   Program  EJECT  1 
used  constant  polytropic  efficiencies  determined  from  test 
data  for  the  secondary  air  and  diffuser  and  an  assumed 
drive  nozzle  polytropic  efficiency  of  0.975,  EJECT  2  used 
Eckert's  polytropic  efficiencies,  and  EJECT  3  used  EJECT  1 
efficiencies,  except  that  the  diffuser  efficiency  was 
varied  with  drive  nozzle  total  pressure  from  0.90  to  0.93. 
The  results  of  the  predictions  for  drive  nozzle  I  and  a 
comparison  with  the  test  results  of  run  18  are  shown  in 
Fig.  31  for  4  lbm/sec  secondary  flow,  and  in  Fig.  32  for 
2  lbm/sec.   Appropriate  changes  were  made  in  the  program 
to  reflect  the  new  dimensions  of  drive  nozzle  II.   The 
prediction  compared  to  test  results  is  shown  in  Fig.  32. 

The  results  show  that  Eckert's  assumed  polytropic 
efficiencies  greatly  underestimate  the  exhauster  performance. 
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At  high  drive  nozzle  pressure  ratios,  the  use  of  the  theo- 
retical prediction  program  EJECT  I  with  experimentally- 
determined  exponents  very  closely  establishes  the  actual 
performance  of  the  exhauster.   The  accuracy  of  the  predic- 
tion decreases  as  drive  nozzle  pressure  ratio  decreases. 
At  the  low  pressure  ratios  the  exhauster  performs  much 
better  than  predicted. 

Much  of  the  difference  noted  above  can  be  attributed 
to  the  physical  formation  of  an  oblique  shock  within  the 
drive  nozzle  for  off -design  operation.   The  classical 
approach  to  describing  nozzle  operation  notes  a  third 
critical  pressure  which  is  the  exit  static  pressure  that 
provides  design  operation.   As  the  receiver  static  pres- 
sure is  increased,  an  oblique  shock  forms,  originating 
at  the  exit  plane  edge.   Further  increase  in  static  pressure 
increases  the  angle  of  the  oblique  shock  until  a  normal 
shock  forms  in  the  exit  plane  for  the  second  critical  pres- 
sure.  Still  further  increase  in  pressure  causes  the  normal 

shock  to  move  from  the  exit  toward  the  nozzle  throat.   This 
idealized  flow  development  does  not  occur  in  actuality 
because  of  friction  and  three-dimensional  flow  effects. 
The  latter,  in  particular,  produce;  internal  oblique  shocks 
which  have  the  effect  of  extending  the  shock  action  over  an 

extended  axial  distance  in  the  nozzle.   Tests  conducted  by 

31 
Stodola   show  internal  shock  effects  over  the  axial  distance 

of  a  nozzle. 

31Hall,  N.  A.,  Thermodynamics  of  Fluid  Flow  (Englewood 
Cliffs,  N.  J.:  Prentice  Hall,  Inc.,  1951),  Fig.  9.5,  p.  136. 
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By  removing  the  mixing  pipe  and  operating  as  a  free 
jet,  the  presence  of  an  internal  oblique  shock  in  the  drive 
nozzle  could  be  made  visible  because  of  condensations  of 
the  atmospheric  humidity  in  the  primary  stream.   Figure  33 
is  a  photo  of  the  free  jet  operation.   The  oblique  shocks 
originated  inside  the  nozzle  and  extended  outside  the 
nozzle  to  their  intersection  with  a  near  normal  circular 
shock  in  the  center  area  of  the  nozzle.   A  sketch  of  the 
observed  pattern  is  shown  in  Fig.  34.   With  increased  drive 
nozzle  pressure,  the  entire  shock  structure  moved  towards 
the  nozzle  exit. 

To  determine  the  displacement  of  this  shock  structure 
in  the  drive  nozzle  under  exhauster  test  conditions,  six 
static  pressure  taps  were  arranged  in  drive  nozzle  I  in 
positions  shown  by  Dl  through  D6  in  Fig.  35.   At  a  secon- 
dary flow  of  4  lbm/sec  the  drive  nozzle  total  pressure  was 
set  at  values  of  PTN/Patm  of  2.65,  2.45,  2.18,  1.90,  and 
1.63  to  duplicate  conditions  of  run  18.   The  measured 
static  pressures  at  the  nozzle  walls  are  shown  in  Fig.  36. 
They  show  that  the  internal  shock  moves  toward  the  nozzle 
exit  if  the  pressure  ratio  is  increased.   It  should  be 
noted  that  for  each  condition  the  pressure  at  the  discharge 
of  the  secondary  nozzle  was  below  that  required  for  oblique 
shock  formation  at  the  exit  plane  for  isentropic  flow 
theory;  that  is,  the  pressure  was  always  less  than  the  2 
critical  pressure. 
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During  the  free  jet  operation  it  was  observed  also  that 
the  primary  jet  did  not  completely  fill  the  exit  plane  when 
oblique  shocks  occurred  inside  the  nozzle.   There  existed  a 
slowed  down  region  of  flow  at  the  outer  perimeter  of  the 
nozzle  discharge  that  had  a  radial  thickness  of  0.5  to  0.75 
inch.   This  phenomenon  can  be  noticed  in  Fig.  34  and  must 
be  responsible  for  the  apparent  thickness  of  the  mixing  zone 
at  the  nozzle  exit  which  is  illustrated  in  Fig.  24. 

The  order  of  magnitude  of  expected  boundary  layer 
thickness  at  the  nozzle  exit  was  determined  with  Eq.  22 

which  can  be  used  to  evaluate  the  boundary  layer  thickness 

32 
on  a  flat  plate  in  turbulent  flow.     The  Reynolds  number 

applied  is  based  on  distance  from  the  nozzle  throat  to  its 

exit. 


J     =    0,31  (Z£)  (22) 


&1 


6    -  boundary  layer  thickness 

W  *  kinematic  viscosity  (ft/sec2) 

I    -   distance  from  throat  to  exit  (ft) 

V  ■  velocity  (ft/sec) 

Based  on  average  conditions 

£  £  0.07  inch 
This  boundary  layer  thickness  is  much  smaller  than  the  ob- 
served distance  of  about  0.5  to  0.75  inch.   The  increased 


32Schlichting/  H. ,  Boundary  Layer  Theory  (New  York: 
McGr aw  Hill  Book  Co.  Inc.,  1955),  Eq.  2.9,  p.  38. 

54 


boundary  layer  thickness  must  be  caused  by  flow  separations 
from  the  nozzle  walls.   This  separation  is  the  result  of 
the  internal  oblique  shock  structure  that  exists  in  the 
drive  nozzle.   The  nearer  the  throat  the  separation  occurs, 
the  less  of  the  diverging  portion  of  the  nozzle  is  used  to 
increase  Mach  number  of  the  flow;  that  is,  full  expansion 
is  not  realized.   The  effect  is  readily  apparent  by  ob- 
serving average  nozzle  Mach  numbers  obtained  at  station  Ml 
for  pTN/Patm  of  2.65,  2.18,  and  1.63.   Test  data  showed 
Mach  numbers  of  2.08,  1.18,  and  0. 88  respectively.   Figure  37 
is  a  sketch  of  the  shock  conditions  that  appear  to  exist. 

The  effect  of  the  internal  oblique  shock  and  resultant 
primary  flow  separation  is  to  raise  the  static  pressure  of 
a  portion  of  the  nozzle  exit  area  to  the  pressure  of  the 
static  pressure  at  the  discharge  of  the  secondary  nozzle. 
The  momentum  theorem  expressed  by  Eq.  17  uses  in  the 
actual  prediction  program 

PSUM  A9   =  P_„  A,  +  P1m  Alm 
*     IN  1    IT  it 

A2   =  mixing  pipe  area 

A,   =  secondary  flow  annulus  area 

A>l  =  nozzle  exit  area 

P,„  =  nozzle  exit  plane  static  pressure 

P,   =  secondary  air  static  pressure. 

In  the  program  the  value  of  PSUM  is  increased  by  itera- 
tion until  the  conditions  at  the  diffuser  exit  match  the 
atmospheric  pressure.   In  the  prediction  program  the 
assumption  is  made  that  P,N  and  P^T  are  constant, 
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respectively,  over  A!  and  *»-  .   In  reality  P^  occurs 
on  an  area  smaller  than  A1#  and  P1T  on  an  area  greater 
than  A^T/  for  conditions  where  the  internal  oblique  shock 
exists.   The  result  is  that  for  a  given  value  of  PSUM  A2, 
the  pressure  PlT  may  be  smaller  than  predicted  to  satisfy 
the  equality.   Thus,  the  hood  total  pressure  would  also  be 
lower.   This  condition  should  be  taken  care  of  in  the  pro- 
gram for  a  more  correct  performance  evaluation  at  low 
pressure  ratios. 

11.   Conclusions  and  Recommendations. 

The  performance  analysis  of  the  exhauster  system  based 
on  polytropic  efficiencies  and  conservation  of  momentum 
very  closely  predicts  the  exhauster  performance  with  re- 
gards to  hood  total  pressure  at  high  drive  nozzle  pressure 
ratios.   The  prediction  under -estimates  exhauster  perform- 
ance at  the  low  pressure  ratios  because  of  the  presence  of 
internal  oblique  shocks  in  the  drive  nozzle.  Accuracy  ob- 
tained by  the  prediction  program  was  the  same  for  drive 
nozzle  IX  which  had  dimensions  different  from  nozzle  I. 

The  hood  total  pressure  is  significantly  effected  by 
the  presence  of  a  traverse  probe  in  the  mixing  pipe.  The 
amount  of  hood  pressure  rise  is  a  function  of  probe  posi- 
tion. Therefore,  traverse  information  presented  cannot  be 
considered  as  having  been  obtained  for  precisely  the  same 
operating  conditions. 
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Traverse  profiles  near  the  drive  nozzle  exhibit  steep 
gradients  of  pressure,  temperature,  and  velocity.   Velocity 
profiles  obtained  appeared  irregular.   The  lower  centerline 
velocity  is  probably  due  to  the  shock  structure  that  occurs 
in  the  nozzle.   It  appears  that  a  sufficient  condition  to 
determine  the  start  of  the  main  mixing  region  in  a  restric- 
ted jet  is  to  determine  where  the  velocity  profile  may  be 
non-diraensionalized  to  a  single  profile.   Results  from  flow 
traverses  showed  that  the  initial  mixing  zone  expands 
radially  at  a  faster  rate  than  predicted  by  Victorin's 
incompressible  approximation. 

Secondary  flow  polytropic  efficiencies  are  independent 
of  drive  nozzle  pressure  but  dependent  on  secondary  flow 
rate.   Diffuser  polytropic  efficiency  is  independent  of 
flow  rate  and  slightly  dependent  on  drive  nozzle  pressure. 
The  diffuser  polytropic  efficiency  is  about  0.90  which  is 
considerably  higher  than  the  previous  estimate  of  0.70. 

Little  test  data  are  available  on  the  behavior  of  com- 
pressible free  jets.   The  facility  used  for  the  ejector 
tests  could  be  adapted  to  obtain  test  data  for  experiments 
in  this  area.   In  addition,  further  insight  into  the  shock 
pattern  in  the  drive  nozzle  and  the  mixing  pipe  might  be 
gained  if  an  adaptation  for  Schlieren  photography  were  de- 
signed.  A  study  of  the  effects  of  the  frictional  and 

three-dimensional  effects  on  shock  formation  in  a  circular 
nozzle  is  considered  the  next  step  necessary  to  improve  the 

exhauster  prediction  analysis  so  that  it  produces  more 
accurate  results  for  the  entire  operating  range. 
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ANNULUS 
STATIC  PRESSURE 
(PSIA) 

6.47 

6.42 

6.37 

6.32** 

6.27 

6.22 

6.17 


POLYTROPIC  EXPONENT 

Secondary  Flow 
(nd) 

1.47 
1.39 
1.34 
1.30 
1.27 
1.24 
1.22 


** 


Actual  Test  Reading 


VARIATION  OF   n 


WITH  ASSUMED  ANNULUS  STATIC  PRESSURE 


WT  =  1.92  Ibm/sec     TTD  =  547.5°      PTD  =6.8  PSIA 


TABLE  II. 
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PTN/patra  nc 

2.65  0.90 

2.45  0.91 

2.18  0.92 

1.90  0.93 

1.63  0.94 


Variation  of  Diffuser  Polytropic  Efficiency 
With  Drive  Nozzle  Pressure  Ratio 


TABLE  IV. 
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Hood  Total 
Pressure 

*Patm  ~  phood' 


Traverse 
Probe  Location 


-6.9  in.  Hg 

-7.0 

-7.2 

-7.5 

-7.85 

-8.15 

-8.55 


Ml 
M2 
M3 
M4 
M5 
M6 
No  Probe 


Hood  Total  Pressure 
versus 
Traverse  Probe  Location 


PTN  =  32  psia 


w.  =  4  lbm/sec 


EJECTOR  TESTS 


TABLE  V. 
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Turbine  Inlet 
Manual  Valve 


Exhauster , 
Plenum 


Diffuser 


ixing    Pipe 


C59"         Rotor    Support 
Bearing 

Shaft  Seal 


Tank  2 


Dynamometer 

/     Dynamometer 
Support  Bearing 


3     M  Exhauster 

lozz,e        /Inlet  Valve 


PIPING     INSTALLATION 
TRANSONIC   TURBINE   TEST  RIG 

FIGURE   I 
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5  teeth  as  on  right  side 


Shaft 


SHAFT    LABYRINTH  SEAL  CONFIGURATION 
TRANSONIC  TURBINE  TEST   RIG 

FIGURE  4 
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Numbered  £  ' s  show  pressure  tap 
locations;  }  taps  equally  spaced 
each  location 


DRIVE    NOZZLE  I  DESIGN 

EJECTOR  TEST  SERIES 

FIGURE  5 
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<t 


6  8   ~ '2  CLASS  2    THREADS 


PITCH  DIA 


6.0709 


^    6  0621 


MAJOR  DIA. 


6  1260 
6  1238 


600 

THO  RQ. 


4  392 


DRIVE   NOZZLE  II    DESIGN 
EJECTOR  TEST    SERIES 
FIGURE  6 
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TRAVERSE     PROBE    INSTALLATION 

EJECTOR  TEST    SERIES 

FIGURE  10 
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APPENDIX  A 

Exhauster  Performance  Analysis 

An  analysis  was  made  of  the  exhauster  performance  to 
predict  hood  total  pressure  which  essentially  parallels  the 
preliminary  design  calculations  of  Vavra.   However,  the 
previous  work  assumed  that  the  nozzle  operates  at  design 
conditions  and  that  the  static  pressure  at  the  entrance  to 
the  mixing  section  is  constant  with  radius.   The  analysis 
explained  in  the  following  permits  a  discontinuity  in  sta- 
tic pressure  between  the  nozzle  exit  (Area  Al)  and  the 
secondary  flow  annulus  (Area  AIT)  at  the  start  of  the  mix- 
ing pipe. 

In  order  to  simplify  the  equations,  polytropic  pro- 
cess relations  as  well  as  non-dimensional  flow  functions 
were  used.   These  were  taken  from  Vavra. 

For  an  expansion  process 


/up™  h         J-e      V*-1  LUJ       \fTI 


ml 


(Al) 


where 


Y 


(A2) 
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n  =  polytropic  exponent 

Y?p  =  polytropic  efficiency        J 

*  =  ratio  of  specific  heats 

w  =  flow  rate  (lbra/sec) 

Tt  =  inlet  stagnation  temperature  (°R) 

•"O 

2 

A2   =  Flow  area  at  station  of  interest  (in.  ) 


P.   =  Inlet  stagnation  pressure  (psia) 
ro 

P2  -  static  pressure  at  station  of  interest 
(psia) 

Tb / 

f  sec^  ' 

Subscript  1  refers  to  inlet  conditions 


Subscript  2  refers  to  station  of  interest 
For  a  compression  process 


•^n? 


w 


fJA; 


t-i-wmrwi 


(A3) 


'2. 

where 

p-^  =  static  pressure  at  1 

Y)       -  *  ^P (A4) 

The  flow  through  the  nozzle  will  be  choked  when  it  is 
operating  at  supersonic  conditions.   Therefore   xe   (or 
xe   )  must  be  a  maximum  value.   Differentiating   }     and 
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setting  the  result  equal  to  zero  yield  the  critical  pres 
sure  ratio  for  choked  conditions 


T'  'cnTlC<xl  x 


Substituting  this  value  into  Eq.  Al,  the  expression  for  the 
flow  function  of  a  choked  nozzle  simplifies  to 


§ 


j 


ZX    /n-1 


(A5) 


chokeJL       \Y)±1    I  W-1    1^  +  1 

Nomle 


It  was  assumed  that  the  drive  nozzle  total  temperature 
and  total  pressure  were  known / as  well  as  the  secondary  flow 
rate  and  total  temperature.   Polytropic  efficiencies  for 
the  drive  nozzle  are  assumed  or  determined  from  testing  and 
held  constant  over  the  operating  range  of  the  exhauster. 
Notation  in  the  development  is  consistent  with  Fig.  3. 

Drive  Nozzle 

For  a  choked  nozzle, 

i 


£ 


2JL  /JVJ- \  (A6) 


\"V  /      y*-i  \r)r  +  \ 


The  flow  functions  from  the  nozzle  entrance  to  the  throat 
(  ?n  ),  and  from  the  nozzle  entrance  to  the  exit  (  $e)<  )  are, 
respectively, 
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I 


*N       Vtn 


n 


AT   PT 


TN 


IF 


5 


wn"yttn  *  inr 

A,    PTn      If  3 


(A7) 


(A8) 


There  are 


Xy>  AT  IwM| 


(A9) 


and 


/V 


*W      =        i     /I 


(A10) 


Expressing  the  exit  flow  function  in  terms  of  the  pressure 
ratio  P^/P^ 


h>m(tr-(m 


(All) 


With  ?ex  given  by  Eq.  A10,  and  with  PTN  known,  PIN  can  be 
solved  by  iteration.   Then 


(A12) 


V,n     s  V^lJ<:P(VTN-TIN) 

w      =   i     p   a  V"^1 

N      xex  'tn  H>  Vft-r 


(A13) 


(A14) 
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Secondary  Flow 

Secondary  air  is  that  part  of  the  high  pressure  air 
provided  by  the  Allis-Chalmers  compressor  that  passes 
through  the  test  turbine  to  the  hood.   For  the  expansion 
process  from  the  hood  to  the  annul us  area 


A.t  Ptd 


If*-  luJ  -tar. 


(A15) 


With  the  secondary  flow  rate  and  hood  total  temperatures 
known,  a  value  of  PIT  was  assumed,  and  PTD  is  calculated 
by  iteration  and  matching.   For  the  assumed  PIT 


T,    =  TTD 


Rr 

Pro 

v,T  -  Yl^cTnvo  -  t1t) 


(A16) 


(A17) 


Mixing  Pipe 

The  momentum  equation  was  used  for  the  mixing  pipe  to 
obtain  conditions  at  the  diffuser  entrance.   The  flow  rate 
through  the  drive  nozzle  is  known  from  preceding  calcula- 
tions and  values  must  be  assumed  for  the  secondary  flow 
conditions.   Applying  the  momentum  equation  to  the  control 
volume  between  the  nozzle  exit  and  annulus  plane  (sta- 
tion 1),  and  the  entrance  to  the  diffuser  (station  2) / 
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/wn+wt)v   _  *W    _  v^r  w-    -  p  ft     +  p    ft   -pa    _C  (A18) 

\    3      /  *■     IP  T^    T      T  *      f 

Equation  A18  holds  for  uniform  velocity  profiles  and  uniform 
static  pressure  distribution  over  each  area. 

The  friction  force  Ff  may  be  expressed  as  a  pressure 
loss 

Ff  =  T  DZ  *V   =   *?  A*-  <A19> 


where 


"=4fV'  (A20) 


L  ■  length  of  duct 
D2  ■  duct  diameter 

V  ■  average  velocity 

JO    ■  average  density 

f  =  friction  factor 

The  friction  factor  was  determined  to  be  0.017, in  the  man- 

33 
ner  outlined  by  Naviaux.    The  friction  force  was  figured 

with  one  iteration  only.   For  the  initial  solution  of 
Eq.  A18,  Ff  was  considered  zero.   After  obtaining  the  con- 
ditions at  the  diffuser  entrance,  average  values  of  V  and 
~p     for  the  mixing  pipe  were  calculated, and  a  A  P   found 
using  Eq.  A20.  AP  A2  was  then  subtracted  from  ^1^1^  + 
P1NA1N   to  obtain  a  new  solution  Of  Eq.  A18. 


33Naviaux,  op.  cit. ,  pp.  25  -  26. 
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Unknown  quantities  in  Eq.  A18  are   V2   and  p2.   The 
pressure  p2  can  be  expressed  by  V2   to  obtain  one  inde- 
pendent variable/  if  the  total  temperature  after  mixing  is 
taken  as  the  mass  weighted  average  of 


T 


*r  Tin    +  ^n  TrN 


WN  +  *V 


(A21) 


This  condition  assumes  the  existence  of  an  adiabatic  pro- 
cess.  With 


wXT  +  wm  =  w 
N    T 


(A22) 


and  the  continuity  equation 


W    = 


*  K\ 


(A23) 


v*  = 


vRT2 


w 


A«1 


w 


v  R 


«**>* 


Tr^    ~ 


(A24) 


Therefore 


■p     __    w  R 


M*     L 


T, 


T2 


-A'   ] 


(A25) 


Substituting   Eq.    A25    into   Eq.    A18 


(A26) 
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Multiply  Eq.  A26  by  V2  and  rearrange 


13     z^Tcp 


Y,    - 


H+fX^A^A.k  +  WRTT^  =  0 


(A27) 


Let 


Bl   = 


?''- 


aJc 


(A28) 


B2    =- 


Stew    -fJttV  +  P  A   ■+  P  A 

rviN    *  T|T*nr«iT^  'i^i 


3    ''*     3 


(A29) 


B3   = 


vvRTJ.; 


(A30) 


Bl   V2  +   B2   V2  +   B3   =   0 


(A31) 


Thus 


V 


-P*-'VBi'-Hg,BJ  _    0Z 


2  6, 


2  6, 


IT1 1 


(A3  2) 


and 


T     =   T       -  J4l 


p     =    p   (III 


(A3  3) 


(A34) 


The  negative  sign  before  the  square  root  yields  the 
subsonic  solution.   Supersonic  conditions  at  the  diffuser 
are  outside  to  operating  range  of  the  TTTR  installation. 


110 


Finally,  for  the  compression  in  the  subsonic  diffuser 


tf  w  - 1  -wmw 


(A3  5) 


P3   is  solved  by  iteration  and  compared  to  atmospheric 
pressure.   If  it  does  not  match  for  the  chosen  values  of 
secondary  flow  conditions,  a  new  PIT  is  selected  and  the 
process  repeated  until  p^   equals  P^tm*   With  P3 
obtained 


«c-' 


T3   =    T2l(f^< 
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(A3  6) 


V3   =  Y*lJe,(Tn-Tai) 

PT3   =  f3    (  xT  ]   '  «A38) 


APPENDIX  B 

Exhauster  Performance  Prediction  Program 

The  analysis  of  exhauster  performance  presented  in 
detail  in  Appendix  A  was  programed  in  Fortran  IV  for  cal- 
culation on  the  IBM  360  computer.   Input  data  include  the 
assumed  polytropic  efficiencies,  drive  nozzle  total  pres- 
sure, drive  nozzle  total  temperature,  secondary  air  total 
temperature,  secondary  flow  rate  and  atmospheric  pressure. 
The  principal  quantity  to  be  determined  is  the  predicted 
hood  total  pressure.   The  printed  output  does  include 
final  values  of  temperature,  velocity  and  pressure  calcu- 
lated during  the  computations. 

The  program  uses  an  executive  MAIN  routine  to  arrange 
the  sequence  of  calculations/  but  all  calculations  are 
performed  in  subroutines.   The  name  and  function  of  each 
subroutine  is: 

CSTANT;  a  grouping  of  all  the  basic  quantities 
used  in  the  entire  program. 

CHOPHI;  a  calculation  of  the  non-dimensional  flow 
function  for  a  choked  nozzle. 

NOZZLE?  computes  conditions  at  the  nozzle  exit. 

SECFLO;  computes  conditions  at  the  secondary  flow 
annulus  area  at  the  entrance  to  the  mixing  pipe. 

MMENTMj  Using  conditions  determined  by  NOZZLE  and 
SECFLO,  conditions  at  the  diffuser  entrance  are  calculated 
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based  on  conservation  of  momentum.   One  iteration  is  in- 
cluded to  determine  the  wall  frictional  shear  force. 

SBCOMP;  computes  the  compression  process  in  the 
diffuser  to  obtain  conditions  at  the  diffuser  exit. 

Input  information  required  to  operate  the  prediction 

program  is: 

Card  1; 

ETAD   -  Secondary  flow  polytropic  efficiency 

ETAN  -  Nozzle  polytropic  efficiency,  entrance  to 
throat 

ETA1   -  Nozzle  polytropic  efficiency,  entrance  to 
exit 

NSETS  -  Number  of  sets  of  data  to  be  computed 

Card  2; 

PTN   -  Drive  nozzle  total  pressure  (psia) 

TTN   -  Drive  nozzle  total  temperature  (°R) 

TTD   -  Hood  (secondary  air)  total  temperature  (°R) 

WT    -  Secondary  flow  rate 

PA    -  Atmospheric  pressure 

Program  EJECTOR  with  a  sample  output  sheet  is  given  on  the 

following  pages. 
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APPENDIX  C 
Data  Reduction  Programs 

Data  obtained  from  laboratory  tests  to  compute  secon- 
dary flow  rate,  velocity  profiles,  and  polytropic  exponents 
were  reduced  by  the  computer  programs  included  in  this 
section.   Input  requirements  are  listed  on  comment  cards  at 
the  start  of  each  program. 

Program  SECFLO  computes  the  mass  of  air  that  passes 
from  the  hood  to  the  mixing  pipe.   The  program  basically 
calculates  Eq.  1  of  the  main  text.   In  addition,  the  small 
quantity  of  air  that  enters  the  hood  through  the  shaft 
labyrinth  is  approximated  by  Eq.  2.   Naviaux's  flow  nozzle 
coefficient  is  determined  as  a  function  of  Reynolds  number 
by  a  series  of  "IF"  statements.   The  values  used  range  from 
1.00  -  1.05  in  increments  of  .01. 

Program  VELOCITY  computes  Mach  number  and  Velocity 
from  traverse  probe  data.   Input  pressures  are  tested  for 
supersonic  or  subsonic  conditions.   Either  Eq.  4  or  Eq.  5 
is  used  depending  on  the  condition  that  exists. 

Program  POLY  computes  by  iteration  the  secondary  flow 
polytropic  exponent  from  test  data.  Eq.  19  with  all  quan- 
tities known  but  nd  was  used  in  the  calculation. 
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APPENDIX  D 
Restricted  Jet  Operation 

In  order  that  the  reader  might  better  understand  the 
general  terras  used  in  this  thesis  to  describe  jet  operation 
and  associated  mixing  phenomena,  the  following  general 
description  of  jet  operation  is  presented. 

Figure  Dl  of  this  appendix  pictures  the  physical  sit- 
uation to  be  discussed.   Essentially,  two  co-flowing 

33 
streams  enter  a  mixing  pipe.   Abraraovich   states  that  with 

this  situation,  the  velocity  components  along  the  y-axis 
(the  radial  components  for  the  exhauster)  prove  so  small 
that  they  can  be  disregarded  in  engineering  problems  in- 
volving jet  theory. 

The  two  co-flowing  streams  enter  the  mixing  pipe  with 
their  fluid  particles  traveling  at  different  velocities. 
The  viscous  forces  acting  between  the  two  flows  cause  a 
slowing  down  of  the  particles  of  the  primary  jet  and  a 
speeding  up  of  the  secondary  flow  particles.   The  area  of 
the  flow  where  the  interaction  between  the  particles  of 
the  two  streams  takes  place  is  called  the  mixing  zone. 
The  mixing  zone  thickens  with  axial  distance,  causing  a 
gradual  "eating  up"  of  the  primary  and  secondary  core 
regions. 

Two  main  mixing  regions  are  designated  in  the  mixing 
pipe.   The  initial  region  has  the  axial  distance  where  the 


33 

Abramovich,  op.  cit. ,  p.  5 
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secondary  core  region  and  the  primary  core  region  still 
exist.   The  main  mixing  region  is  that  part  of  the  pipe 
where  the  mixing  zone  has  touched  the  pipe  walls  and  the 
primary  core  region  no  longer  exists.   The  regions  are 
shown  in  Fig.  Dl.   In  the  literature  use  is  sometimes  made 
of  a  transition  zone  which  is  the  region  that  occurs  if  the 
mixing  zone  contacts  the  pipe  wall  before  the  primary  core 
is  "eaten  up"  or  vice  versa.   This  thesis  uses  the  simpli- 
fied model  in  which  it  is  assumed  that  the  length  of  the 
transitional  region  is  equal  to  zero.   In  the  main  mixing 
region,  the  jet  becomes  similar  in  appearance  to  a  flow 
of  fluid  from  a  source  of  infinitely  small  thickness  (in 
an  axially  symmetric  case, the  source  is  a  point). 
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